RECENT EXTENSIONS TO THE FREE- VORTEX- SHEET THEORY 
FOR EXPANDED CONVERGENCE CAPABILITY 


James M. Luckring 
NASA Langley Research Center 
Hampton, Virginia 

Keith D. Hoffler 
Vigyan Research Associates, Inc. 
Hampton, Virginia 

Arthur C. Grantz 
Northrop Corporation 
Hawthorne, California 


SUMMARY 

A new version of the free- vortex- sheet formulation is presented which has 
greatly improved convergence characteristics for a broad range of geometries. The 
enhanced convergence properties were achieved largely with extended modeling capa- 
bilities of the leading-edge vortex and the near-field trailing wake. Results from 
the new code, designated FVS-1, are presented for a variety of configurations and 
flow conditions with emphasis on vortex flap appl ications. 


INTRODUCTION 

The design constraints for high-speed aircraft with efficient supersonic 
performance capability often result in highly swept wings which are conducive to 
separation-induced leading-edge vortex flows at moderate-to-high angles of attack. 
Although this type of separation can cause increased drag through the loss of 
leading-edge thrust, the vortex-lift characteristics resulting from this highly 
stable flow can be used to improve takeoff and landing capabilities, provide 
maneuver lift increments at subsonic and transonic speeds, and, most recently, 
reduce drag at moderate to high lift coefficients through the use of simple "vortex 
flaps." As for attached flow, however, the full exploitation of this vortex flow 
will require improved theoretical analysis and design techniques capable of pro- 
viding reasonable estimates of the three-dimensional surface pressure distributions 
at a practical expenditure of human and computer resources. 

However, the computation of these vortex flow effects has proven to be a 
challenging task. It is only during the last 10 years or so that methods have begun 
to emerge with the stated capability, even for simple three-dimensional geometries 
with sharp leading edges at subsonic speeds. Among these methods, the free-vortex- 
sheet theory has provided the best estimates of the inviscid surface pressure 
distributions for a rather broad class of generic wing shapes. The major drawback 
of this method has proven to be the difficulty often encountered in achieving 
converged results. 



C, A new version of the free- vortex- sheet formulation, designated FVS-1, is 
presented in this paper which has greatly improved convergence properties for a 
broad range of geometries including vortex flaps. The convergence difficulties were 
approached from a fluid-mechanical viewpoint as opposed to that of numerical 
analysis. As a consequence, the basic numerical method for iteratively solving the 
nonlinear vortex flow problem required little change. The enhanced convergence 
properties were achieved largely with extended modeling capabilities of the leading- 
edge vortex and the near- field trailing wakes. Using this approach, considerable 
insight was gained into the underlying fluid mechanics of vortex modeling as related 
to convergence considerations. Results are presented which demonstrate the expanded 
convergence capabilities with emphasis on the salient features modeled as well as 
the aerodynamic consequences of these modeling techniques. 


SYMBOLS 


AR aspect ratio 

b wing span 

Cb static root bending moment coefficient, static root bending moment/ 

(q« Sref ^) 

Cq drag coefficient, drag/Iq,*, S re f) 

Cq 9 o experimental value of drag coefficient at Cl = 0 
aCq drag-due-to-lift coefficient, Cq - Cd,o 

Cl lift coefficient, lift/iq*, S re f) 

C m pitching moment coefficient, pitching moment/fq^ S ref c) 

Cfj normal force coefficient, normal force/lq*, S re f) 

Cp static pressure coefficient, (p - Pco)/q ro 

ACp Cp jU - Cp } i 

c streamwise chord 

c reference chord 

c r root chord 

ct tip chord 

E log (e r 2 ) 

It iteration number 

i nw longitudinal extent of near wake in fraction of vortex-core diameter at 

trailing edge 
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Mach number 

number of update iterations 
unit normal vector 

static pressure; also maximum camber height in fraction local semi span 

p 

freestream reference dynamic pressure, (1/2) U„ 

individual panel residual error 

reference area 

local semi span 

thickness 

freestream reference velocity 
perturbation velocity components 
body-axis Cartesian coordinates 
angle of attack, degrees 
angle of attack for smooth onflow 


Aa a ~ ago 

p angle of sideslip, degrees; also, Prandtl-Glauert factor, (1 - M* 2 ) 1 / 2 

6 deflection angle, degrees, positive downward 

5 ^ vortex flap deflection angle normal to hingeline, degrees, positive 

downward 

A6f increment in 6^ 

6 t side-edge rake angle, degrees, positive inboard 

A sweep angle, degrees 

\ taper ratio, ct/c r 

d freestream reference density 

y oo 

E summation 

t vorticity vector 

. vector dot product 

)| vector magnitude 
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Sub sc ri pts 
f fl ap 

hi hingeline 

1 lower surface 

le leading edge 

te trailing edge 

u upper surface 

1 first order in perturbation quantities; also, inboard 

2 second order in perturbation quantities; also, outboard 

oo free stream reference conditions 

Abbreviations 

Arg argument 

FVS free vortex sheet 


CONVERGED FVS PARAMETER SPACE 

A study was initially conducted to document the converged parameter space of 
the free-vortex-sheet method as demonstrated in the open literature. The results 
are detailed in Table 1 and partially summarized in figure 1. In Table 1, the 
results are organized by planform type, first for flat wings and then for cambered 
wings; within one class of planform and camber they are listed alphabetically by 
first author. It should be noted that most entries correspond to a range of 
solutions obtained for a systematic variation of some salient parameter such as 
angle of attack, Mach number, sweep angle, flap deflection angle, etc. Although a 

number of significant aspects of this solution space could be called out (not to 
mention details of the solutions themselves), the main point to be observed is the 
extent of geometric variation and freestream conditions over which the method has 
been successfully applied. 

Some general features of free-vortex- sheet solutions are highlighted in figure 
2. By virtue of explicitly modeling the leading-edge vortex, a variety of "real 
flow" effects are included in the FVS solutions which, to date, have not been 
calculable by simpler methods. The chief feature of these solutions is an accurate 
estimation of the three-dimensional, inviscid pressure field. By way of integration 
this also results in accurate force and moment estimates. In addition, high angle- 
of-attack vortex lift loss effects due to longitudinal vortex curvature near the 
trailing edge as well as vortex crowding are implicitly represented. Additional 
discussion of these effects has been given by Polhamus (1983, 1985). 

Although not covered in this paper, the FVS method has recently been extended 
to account for two primary effects of viscosity and distributed vortici ty , vortex 
breakdown and secondary vortex separation. (The method is still limited to flows 


88 

C ' ^ 



with a known primary separation line such as occurs with a sharp leading edge.) 

By imbedding a parabolized Navier-Stokes representation of an axi symmetric vortex 
core, Luckring (1985) showed that the condition of incipient vortex breakdown at the 
trailing edge of delta, arrow, and diamond wings could be predicted with a critical 
helix angle concept. Wai et al . (1985) demonstrated the first estimations of 
secondary vortex separation with the FVS method by modeling the upper surface 
boundary layer flow and interacting the viscous flow effects with the outer inviscid 
flow. Although continuing development of the vortex core and secondary vortex 
effects will be required, the initial findings of these studies are quite promising. 
The current status of secondary separation modeling by boundary layer techniques 
will be addressed in this conference by Blom et al . (1985) as well as by 
Woodson and DeJarnette (1985). 


METHOD EXTENSIONS 

From the experience gained in utilizing the free-vortex- sheet method, four 
aspects of the formulation were chosen for modification. The first three aspects 
regard flow modeling considerations; they are the near-wake model, the restart 
capability, and the manipulation of the starting vortex sheet geometry. The fourth 
aspect regards the quasi-Newton scheme employed to solve the nonlinear system of 
governing equations. 


Near-Wake Considerations 

The primary attributes of the near-wake flow field are summarized in figure 3 
for the conditions of planar and nonplanar vorticity. In either case, the trailing- 
edge Kutta condition results in a near-wake flow which is tangent to the wing at the 
trailing edge. For planar vorticity the linear form of the pressure coefficient is 
appropriate and, when taken in conjunction with the no-load wake boundary condition, 
results in a wake flow where the magnitude of the vorticity vector is free, but its 
argument is fixed at the streamwise direction. 

As is true for slender-body theory, nonplanar vorticity necessitates the 
inclusion of nonlinear terms in the equation for the pressure coefficient. Although 
these terms are second order in form, their contribution is of first order, particu- 
larly with regard to the si dewash (v). Employing this more exact form of the 
pressure coefficient in conjunction with the no-load wake boundary condition yields 
a more realistic wake flow where both the magnitude and the argument of the 
vorticity vector are now free quantities. These effects allow the wake vorticity to 
skew laterally and significantly affect the satisfaction of the trail ing-edge Kutta 
condition as was shown by Luckring et al . (1982), for leading-edge vortex flows. 
However, it must be emphasized that the driving mechanism of this wake effect is the 
presence of nonplanar vorticity, only one example of which is the leading-edge 
vortex problem. Other examples of conditions where these effects would also be 
important include wings with significantly nonplanar geometries such as wings with 
winglets or vortex flaps, regardless of whether the flow is attached or separated. 

An example of the nonplanar effect is presented in figure 4 in the form of 
vorticity contours. Although both calculations utilized the higher order near wake, 
the attached flow case of figure 4a resulted in the wake vorticity orienting itself 

largely in the streamwise direction since the vorticity of this flow is planar. 

Here, it is evident that the conventional approximations made for the simple 
trailing wake are adequate for this flow. In the nonplanar vortex- flow case of 
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figure 4b the lateral skewing of the wake vorticity is significant; modeling this 
flow with a simple trailing wake would constitute a poor approximation. It is of 
interest to note the coalescence of vorticity downstream in the near wake. This is 
indicative of the wake vortex which occurs on slender wings with leading-edge 
separation; it is shed along the trailing edge and has the opposite sense of the 
leading-edge vortex. 

From a convergence viewpoint, it is desirable to make the near wake as short as 
possible while maintaining the beneficial effects of this model. A long near wake 
would not only require more panels (increasing computational expense) but also 
require increasingly complicated wake flow effects to be represented, hampering 
convergence. Fortunately, the results of a near-wake length study (figure 5) show 
very little effect of the length of near wake modeled on the converged force and 
moment properties of a delta wing. In this figure, the near-wake length is 
presented in fractions of the vortex core diameter at the wing trailing edge. Also 
shown are the suction analogy estimates of the normal force coefficient for two 
angles of attack. Apparently the salient feature of the higher order near wake is 
to allow the wing vorticity to pass onto the wake with the proper orientation. As 
has been shown in other studies, details of the wake rollup have little effect on 
the wing flow itself. Additional calculations (not shown) indicated little sensi- 
tivity of the solution to near-wake paneling variations at a fixed near-wake length. 

A comparison of the convergence histories for a wing with no near wake and a 
wing with a short near wake is shown in figure 6. Both solutions required eight 
iterations to drive the residual error below a tolerence of 10" 6 , and the overall 
character of the convergence was unchanged by adding the short near wake. It is 
noteworthy that the consequences of modeling nonplanar vorticity effects with the 
near wake have a larger effect on the solution than does driving the residual error 
down 10 orders of magnitude from unacceptable to fully acceptable levels. 

The results of figure 6 also indicate that the overall force and moment 
coefficients are converged at a residual error level on the order of 10~ 3 , as 
opposed to the default tolerance of 10~ 6 . In figure 7, the effect of the residual 
error level on the spanwise pressure di stributions very near the trailing edge 
(where convergence is generally slowest) also indicate that 10" 3 is acceptable for 
practical applications. Examination of the other calculations of this paper further 
supports this conclusion. 

Computed force and moment properties for the uni t- aspect- ratio delta wing of 
Hummel (1979) are compared with experimental values in figure 8. The wing was 
represented by mean plane doublet networks to account for the asymmetric beveling 
and thickness effects were not modeled. The calculations (with a short near wake) 
show excellent agreement with the experiment including camber effects which, for the 
most part, are due to the trail ing-edge bevel (a trai ling-edge- flap type of effect). 
Solutions were obtained down to relatively low angles of attack (in magnitude) and 
lift coefficients; at ±5 degrees the computed vortex-lift increment was approxi- 
mately ±0.03 from the computed attached flow value. (The attached flow solutions 
are for zero leading-edge suction.) 

Spanwise pressure distributions for 20° angle of attack correlate reasonably 
well with the experimental results (figure 9). Differences on the upper surface in 
peak suction magnitude and in the outboard values are largely due to secondary 
separation, a viscous flow effect. At the aft station shown., thickness effects also 
affect the quality of the correlation. This is the station of maximum thickness, 
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and the thin wing calculations reflect the singular nature of the modeled flow at 
this station. 


Partial Restart 

A quasi-Newton method is used to solve the nonlinear equations for the free- 
vortex- sheet formulation. As with any Newton method, the starting solution must be 
"sufficiently close" to the converged answer if that answer is to be achieved. The 
basic FVS code has a full restart capability which is useful for changes in the far- 
field boundary conditions. Here the complete solution (singularities and vortex 
geometry) from a prior calculation is used as a starting solution for a new 
calculation of the same wing geometry (and paneling) at different freestream 
conditions. This capability has proven to be very useful and cost effective for 
computations over a range of freestream conditions such as angle of attack, Mach 
number, etc., and was used for the results of the Hummel delta wing just discussed. 
For cases where no prior information is available, the default approach was to base 
the initial vortex- sheet geometry on the conical -flow solution of Smith (1966) and 
to calculate the initial singularity distribution from a subset of the full 
governing equations. 

The restart capability of this method has been extended to include a "partial 
restart" which is useful for changes in the near-field boundary conditions such as 
would arise for a change in wing geometry, paneling, etc. The partial restart 
allows for the transfer of an arbitrary vortex sheet geometry from some prior 
calculation to the current configuration of interest and solves for the initial 
singularity distribution by satisfying a subset of the full governing equations. 

This technique has been exploited for many of the vortex- flap solutions summarized 
in Table 1, primarily for achieving solutions over a range of vortex- flap deflection 
angles. However, it was found that the utility of the partial restart diminished as 
the differences in boundary conditions between an available and the desired solution 
increased. For this reason, additional vortex-sheet manipulation capability was 
implemented. 

Vortex Sheet Manipulation 

For many practical applications the geometric properties of the leading-edge 
vortex can be dominated by some salient geometric characteristic of the wing. One 
example is the vortex flap application where the resultant vortex geometries are 
similar in a coordinate system normal to the flap surface. Therefore, a simple 
vortex manipulation capability was implemented which allows the initial vortex 
geometry to be arbitrarily rotated and/or scaled about the wing leading edge. This 
idea is illustrated in figure 10, where the vortex sheet from an available three- 
dimensional solution for a planar delta wing is first transferred (via the partial 
restart) to a wing with a deflected vortex flap, and is then rotated by the flap 
deflection angle. In practice, it is often advantageous to also slightly scale the 
resultant initial vortex geometry as well. The effects of this approach on conver- 
gence are shown in figure 11a. Use of the simple transfer, rotation, and scaling 
readily yielded a converged solution for a large increment in flap deflection; for 
small increments in flap deflection either the transfer alone or in conjunction with 
manipulation readily yielded the converged solution. Without these techniques 
converged results had not been achieved. 

An example of the scaling effects is illustrated in figure 11b. The converged 
solution for a delta wing was transferred to the gothic wing and then scaled to 
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approximate the nonconi cal effects due to the gothic wing pi an form. Once again the 
calculation readily converged; using conical flow information for the starting 
solution, the prior version of the code would not yield a converged result. 

The partial restart and vortex manipulation can be used to obtain solutions for 
more complex geometries and flow conditions. Shown in figure 12 is an example for a 
wing-body configuration with a gothic planform vortex flap deflected at a large 

angle. By organizing the calculations roughly as shown one can obtain intermediate 
solutions which are useful for interpreting the complete configuration solution as 
well as for determining incremental effects of geometric variations. In figure 13 
an example is shown for a similar wing geometry but at sideslip conditions. Drawing 
upon simple-sweep concepts, symmetric calculations were first performed which, by 
the partial restart, provided good starting solutions for the complete 
configuration. 

The approach taken for this case can also be cost effective since the gross 
solution features are established with the two symmetric calculations which 
nominally require one-half the computer resources of a full-span calculation for the 
same paneling and number of iterations. The symmetric precursor calculations could 
have been performed at reduced panel ings, further affecting cost reduction. This 
feature has been exploited for various problems by initiating a solution with a 
coarse paneling and then, using the partial restart, obtaining solutions for finer 
panelings, much as is done with mesh sequencing in finite-difference formulations. 

Modified Quasi-Newton Scheme 

With the default quasi-Newton scheme of the free- vortex- sheet method the 
Jacobian matrix is fully formulated at the starting solution and then only every 
third iteration thereafter. For the intermediate iterations the values of the 
Jacobian are obtained with an approximate update procedure due to Broyden (1973). 

By this approach, there is no need to fully reevaluate the partial derivatives 
comprising the elements of the Jacobian for each iteration. So long as the method 
converges, this technique can reduce the cost of obtaining the solution. As with 
any approximation technique, however, the update method will work best in conditions 
where the elements being approximated vary linearly and, in actuality, undergo 
little change. 

In general, the elements of the Jacobian will undergo greater change early in 
the iteration process, where residual error levels are high, than they will later in 
the calculation as the solution becomes established. Strictly reformulating the 
Jacobian matrix every three iterations can therefore result not only in unnecessary 
numerical difficulties early in a calculation (due to inherent approximations 
affiliated with the update procedure) but also in unnecessary numerical expense 
late in a calculation (due to Jacobian reformulation when an update would be just as 
good). Therefore, the Jacobian reformulation frequency was adjusted to be a 
function of the residual error level. The following approximate schedule was 
determined from numerical experimentation for the number of iterations employing the 
update procedure as a function of residual error level: 


error range N U p 

E > 0.5 1 
0.5 > E > 0.1 2 
0.1 > E 5 
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Several cases of this study were reanalyzed with the above Jacobian update 
schedule. Although no cases have been found which required this technique to obtain 
convergence, calculations with the adjusted Jacobian update schedule evidenced 
smoother convergence properties than did the baseline calculations. The cost of the 
modified calculation was comparable to the baseline calculation primarily because 
the increased Jacobian reformulation frequency at high residual error levels is 
traded off against the reduced frequency at low levels of error. 


A DIRECT APPROACH FOR CAMBERED WINGS 

The techniques described in this report can also be used directly to obtain 
solutions for cases where no prior information is available for use with the partial 
restart. An example is shown in figure 14 for wings with large camber at freestream 
conditions where the leading-edge vortex is small and acting on the forward- facing 
camber surface. The approach, to be subsequently described, will be first applied 
to the circular-arc coni cal -camber wings of Barsby (1974) and then to a vortex flap 
configuration. 

Attached flow calculations were first performed with the free- vortex- sheet code 
to determine the smooth-onflow angle of attack. For these conical wings, this 
condition occurred along the entire leading edge within an angle of attack of about 
1°. The modeling included the short near wake and the calculations incorporated the 
scheduled Jacobian update procedure. Because the angle of attack differs slightly 
from the smooth-onflow angle, the leading-edge vortex will be small and governed 
primarily by properties near the leading edge. To approximate these effects, the 
initial vortex geometry was based upon the conical flow solution of Smith (1966) for 
the incremental angle of attack beyond smooth onflow and was rotated by the 
transverse wing slope at the leading edge. 

Calculations were performed for a unit-aspect-ratio wing, first with a 
nondimensional camber height of 0.4 for freestream conditions 5° above smooth 
onflow. The calculation readily converged, and some properties of the three- 
dimensional solution are shown in figure 15a at the wing trailing edge. The simple 
starting procedure just described provided a reasonable estimate of the converged 
vortex geometry. For comparison purposes, the flat-plate conical-flow solution is 
also shown. This solution emphasizes the significant effects camber has on the 
vortex solution and also indicates probable convergence difficulties if the flat- 
plate solution were to be used as a starting solution for the cambered wing. 

A more extreme case is shown in figure 15b. From the previous calculation 
(p = 0.4) it was apparent that the described starting procedure resulted in an 
oversized vortex, so the starting vortex geometry for the current calculation 
(p = 0.6) was further scaled down by 50 percent. The calculation readily converged 
and resulted in an extremely small vortex. As a consequence, this solution 
exhibited nearly 90-percent suction, based upon the (usual) attached- flow/zero- 
suction and flat-pl ate-opt^ levels. 

As a final application, this direct approach was applied to a vortex flap 
geometry. The resultant convergence histories (figure 16) show the dramatic effects 
of the new formulation. It. should be noted that convergence is displayed at a 

residual error level of 10 ; as was shown earlier the solution is sufficiently 
established at a residual level of 10-3 f or practical purposes. 
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CONCLUDING REMARKS 


A variety of extensions have been implemented to the free- vortex- sheet 
formulation for the purpose of improving the convergence characteristics of the 

method. By extending the modeling capability of the code significant improvements 
in the convergence characteristics were realized without changing the basic 
numerical scheme. Some results of this study are summarized below. 

Only a short near wake is required to capture the significant effects of 
nonplanar vorticity. The convergence of a calculation with a short near wake is 
nearly as fast as the convergence of a calculation with no near wake at all. With 
either wake representation, the solution properties evidenced little variation 
beyond a residual error level of 10** 3 . 

A "partial restart" was added to the formulation. This allows for the use of 
an arbitrary vortex from a prior calculation as a starting solution. In addition, 
simple vortex manipulation features such as rotation and scaling proved to be very 
effective in enhancing convergence. 

With the extensions described in this report, the utility of the FVS 
formulation has been extended to conditions where the vortex is small. This allows 
for three-dimensional calculations to be easily performed at low angles of attack 
and on configurations with highly deflected leading edges. 
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" 

•• 

» 

10-20 


" 

(1975) 

















Brune, 

Delta 

74 

- 

0 

1.15 

0 

Flat 

_ 

_ 

_ 

Thin 

. 

Full 

5-20 

0 

0.8 

et al. 

















(1977) 

















Buter, 

Delta 

74 

- 

0 

1.15 

0 

Flat 

_ 

_ 

_ 

Thin 

_ 

Full 

10 

0 

0 

et al. 

















(1985) 

















Erickson 

Delta 

55 

- 

0 

2.80 

0 

Flat 

_ 

_ 

_ 

Thin 

_ 

Full 

15 

0 

0.6 

(1983) 


65 

" 


1.87 

" 







" 



•• 






" 

" 

" 

'* 



" 

*■ 

Part 



■* 



75 



1.07 



" 



" 


Full 

" 

" 

" 

Gloss , 

Delta 

63.4 

- 

0 

2.00 

0 

Flat 

_ 

_ 

_ 

Thin 

_ 

Full 

20 

0 

0 

et al. 


70 


" 

1.46 

" 



" 


•• 

•• 

•• 

14 

*• 

•• 

(1976) 

















Hoffler 

Delta 

65 

- 

0 

1.87 

0 

Flat 

_ 

_ 


Thin 

_ 

Full 

16 

0 

0 

(1985) 

















Johnson, 

Delta 

70 

- 

0 

1.46 

0 

Flat 

_ 


_ 

Thin 

_ 

Full 

14 

0 

0 

et al. 


74 

" 

" 

1.15 

" 


" 






5-25 

•» 

0.6 

(1976) 

















Johnson, 

Delta 

63.4 

- 

0 

2.00 

0 

Flat 

_ 

_ 

_ 

Thin 

- 

Full 

20 

0 

0 

et al. 


70 



1.46 

" 

" 



’’ 


•* 

*• 

14 

>• 


(1980) 


74 



1.15 

** 

*■ 







15-30 

•• 

•« 



76 

" 

" 

1.00 

” 

" 


” 

" 


*' 


20 

" 




80 



0.71 

" 

" 





” 

" 

15 

5 




86.4 



0.25 

" 




" 

" 



2-25 

0 




63-86 



2-0.25 

" 

" 







20 

" 

" 

Kuhlman 

Delta 

63-86 

- 

0 

2-0.25 

0 

Flat 

_ 

- 

_ 

Thin 

_ 

Full 

15 

0 

0 

(1978) 


63-80 

" 


2-0.71 

*' 








20 

’’ 

" 



82.7 


*' 

0.52 



'* 



” 



10-40 

" 

0.8 






" 



" 






10.3 

" 

0.7 



" 



" 



" 






30 

" 

" 










" 

” 




20.8 

" 

0.5-0. 8 

Lamar 

Delta 

74 

- 

0 

1.15 

0 

Flat 

_ 

_ 

_ 

Thin 

_ 

Full 

15-25 

0 

0 

(1978) 


80 


" 

0.71 





" 



" 

10-25 

" 

'* 

Lamar , 

Delta 

74 

- 

0 

1.15 

0 

Flat 


_ 

_ 

Thin 

_ 

Full 

20 

10 

0.2 

et al. 

















(1979) 

















Luckring , 

Delta 

63 

- 

0 

2.04 

0 

Flat 

_ 


_ 

Thin 

_ 

Full 

20 

0 

0.3 

et al. 


74 

" 

" 

1.15 

" 

" 


" 


" 

" 

" 

10-25 

0,5 

0.2 

(1982) 


" 

" 

" 

" 






*• 


" 

20 

0,15 

» 



82.7 

" 

" 

0.52 

" 


" 

” 

" 


■* 

" 

10 

0 

0.99 









" 




" 

" 

20.8 

0 

0-0.99 

Luckring, 

Delta 

65 

- 

0 

1.87 

0 

Flat 

_ 

_ 

_ 

Thin 

_ 

Full 

15 

0 

0 

et al . 

















(1985b) 

















Polhamus 

Delta 

70 

- 

0 

1.46 

0 

Flat 

- 

_ 

_ 

Thin 

_ 

Full 

10-50 

0 

0 

(1983) 


80 



0.71 

" 




" 

" 

" 

" 

5-50 

'* 

" 

Polhamus 

Delta 

70 

- 

0 

1.46 

0 

Flat 

_ 

_ 

_ 

Thin 

_ 

Full 

10-50 

0 

0 

(1985) 


76 

” 

" 

1.00 

" 

’* 


" 


Thick 

" 


25 

*' 

" 



80 



0.71 



" 



Thin 

" 

" 

5-50 

" 

” 

Reddy 

Delta 

70 

- 

0 

1.46 

0 

Flat 

- 

- 

_ 

Thin 

_ 

Full 

10-30 

0 

0 

(1979) 


76 


” 

1.00 

" 




" 




5-30 

“ 

" 

Reddy 

Delta 

70 

- 

0 

1.46 

0 

Flat 

- 

- 

- 

Thin 

- 

Full 

15-40 

0 

0 


(1981) 



TABLE 1.- CONTINUED 


(b) arrow and diamond wings, uncambered 


Reference 

Planform 

Class 

A i 

a 2 

A te 

AR 

X 

Camber 

Class 


A hl 

h 

Thickness 
Wing Fusl 

Vortex 

Span 

Flow 

a 

p 

M 

Brune, 
et al . 
(1975) 

Arrow 

71.2 

- 

43.6 

2.02 

0 

Flat 

- 

- 

- 

Thin - 

Full 

15.8 

0 

0 

Brune, 
et al. 
(1977) 

Arrow 

74 

~ 

37 

1.47 

0 

Flat 

“ 


' 

Thin 

Full 

5-25 

0 

0.6 

Johnson, 

Arrow 

71.2 

_ 

43 

2.00 

0 

Flat 

_ 

_ 

- 

Thin - 

Full 

11.9 

0 

0.4 

et al. 
(1976) 


74 


37 

1.47 








5-25 


0.6 

Johnson, 
et al . 
(1980) 

Arrow 

71.2 

_ 

43 

2.00 

0 

Flat 




Thin - 

Full 

11.9 

0 

0.4 

Kuhlman 

(1978) 

Arrow 

70 

- 

36.9 

2.00 

0 

Flat 




Thin 

Full 

10-30 

0 

0 

Luckring, 
et al. 
(1982) 

Arrow + 

67 

— 

55 

2.50 

0 

Flat 

' 

' 


Thin - 

Full 

7-14 

0 

0.76 

Polhamus 

(1983) 

Arrow 

70 

- 

36.9 

2.00 

0 

Flat 

- 

- 

- 

Thin - 

Full 

5-50 

0 

0 

Reddy 

Arrow 

70 

- 

36.9 

2.00 

0 

Flat 

- 

- 

- 

Thin - 

Full 

15-40 

0 

0 

(1979) 




28.8 

15.4 

1.82 

1.62 




•• 


.. 


■■ 


■* 

Kuhlman 

(1978) 

Diamond 

70 

- 

-51.3 

1.00 

0 

Flat 

- 

- 

- 

Thin - 

Full 

10-30 

0 

0 

Reddy 

Diamond 

70 

- 

-15.4 

1.32 

0 

Flat 

- 

- 

- 

Thin - 

Full 

15-40 

0 

0 

(1979) 


„ 


-28.8 

-39.5 

-47.7 

-51.3 

1.21 

1.12 

1.04 

1.00 





„ 




.. 

•' 

Reddy 

(1981) 

Diamond 

70 

“ 

-51.3 

1.00 

— 

Flat 

— 

— 


Thin - 

Full 

25-40 

0 

0 


+ B-l Planform 
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TABLE 1 .- CONTINUED 


ORIGINAL PACT- ?8 
OF POOR QUALITY 


(c) cropped wings, uncambered 


Reference 

Planform 

Class 

A 1 

a 2 

Ate 

AR 

Camber 
X Class 


A hl s 2 

Thickness 
Wing Fusl 

Vortex 

Span 

FI ow 
a 

p 

M 

Erickson 

Cropped 

65 

. 

0 

2.80 

0.25 Flat 



Thin - 

Full 

15 

0 

0.6 


(1983) 

Delta 















Erickson 

Cropped 

55 

- 

0 

1.80 

0.2 

Flat 

_ 

_ 


Thin 

Full 

15 

0 

0.4 

(1985) 

Delta 















Luckring, 

Cropped 

63 

- 

0 

0.87 

0.4 

Flat 

- 

_ 

- 

Thin 

Full 

20 

0 

0.3 

et al. 
(1982) 

Delta 















Reddy 

Cropped 

63 

- 

0 

0.87 

0.4 

Flat 

_ 

_ 

- 

Thin 

Full 

15-25 

0 

0.2 

(1981) 

Delta 

80 

" 

" 

0.27 

0.45 







15-35 

” 

0 


Manro 

(1983) 

Cropped 

Arrow 

71.2 - 

43 

1.65 

0.1 

Flat 

" 


" 

Thin 

Thick 

Thin 

Full 

Thick " 

Part 

16 

8 

6 

0 

0.4 

Reddy 

(1981) 

Cropped 

Arrow 

63 

40 

1.07 

0.4 

Flat 

- 

- 

- 

Thin 

Full 

15-25 

0 

0.2 

1 


Reddy 

(1981) 

Cropped 63 
Diamond 

-40 0.74 

0.4 Flat - - - Thin 

Full 15-25 0 

0.2 



(d) 

cranked wings, uncambered 




Reference 

Planform 

Class 

A i 

a 2 

A te 

AR A 

Camber 

Class 

5 i 

A hl 

6 2 

Thickness 
Wing Fusl 

Vortex 

Span 

Flow 

a 

p 

M 


Erickson 

(1983) 

Cranked 

Delta 

72 

62 

0 

1.68 0 

Flat 

- 

- 

- 

Thin 

- 

Full 

15 

0 

0 


Kuhlman 

(1978) 

Cranked 

Delta 

80 

65 

0 

1.60 0 

Flat 

- 

- 


Thin 


Full 

5-30 

0 

0 


Reddy 

(1981) 

Cranked 

Delta 

80 

65 

0 

1.60 0 

0.95 0.2 

Flat 

- 



Thin 

: 

Full 

15-30 

0 

0 


Reddy 

(1981) 

Cranked 

Arrow 

80 

65 

30 

2.08 0 

1.10 0.23 

Flat 

- 

- 

- 

Thin 

- 

Full 

15-25 

15-30 

0 

0 


Reddy 

(1981b) 

Cranked 

Arrow 

70-80 

80 

65 

60-70 

30 

2. 6-2.1 0 
2. 6-1. 6 " 

Flat 




Thin 


Full 

20 

0 

0 
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TABLE 1 


CONTINUED 


(e) miscellaneous wings, uncambered 


Planform Camber Thickness Vortex Flow 

Reference Class K\ Ag A^ e AR A Class A^-j &2 Wing Fusl Span a P M 


Brune, 
et al. 
(1975) 

Gothic 

Arrow 

50 

87 

56.3 

1.60 

0 

Flat 

Erickson 

(1983) 

Gothic 

62.5 

76.6 

0 

1.35 

0 

Flat 

Frink 

(1985) 

Gothic 

60 

74 

0 

1.36 

0 

Flat 

Johnson, 
et al. 
(1980) 

Rectangle 0 

' 

0 

0.2-1, 

.2 1 

Flat 

Reddy 

(1981) 

Raked 

Side 

Edge 

63 

-3-13 

0 

1.07 

0.4 

Flat 


” 

“ 

— 

Thin 

” 

Full 

14.3 

0 

0 

- 

- 

- 

Thin 

- 

Full 

15 

0 

0 

_ 

. 

_ 

Thin 

Thin 

Part 

10-16 

0 

0 

'* 




Thick 

Full 

14 

" 

" 

- 

- 

- 

Thin 

- 

Full 

20 

0 

0.2 




Thin 


Full 

23 

0 

0 


( f ) all planforms, conventional camber 


Planform 

Reference Class A^ A 2 A^ e AR 


Conventional Camber 
X Class 6l Ahl 62 


Thickness Vortex Flow 
Wing Fusl Span a 


M 


Buter, 
et al. 
(1985) 

Delta 

74 


0 

1.15 

0 

Apex Flap 

0-40 

20 

" 

Thin 

Full 

Inbd 

10 0 

0 

Luckring, 
et al. 
(1985b) 

Delta 

76 

“ 

0 

1.00 

0 

T.E. Flap 

-6 

‘ 

Thin - 

Full 

-20-50 0 

0 

Reddy 

Delta 

75 

- 

0 

1.07 

0 

Smooth Onflow, 

- 

- 

Thin - 

Full 

2-30 0 

0 


(1981) Cld « 0.3 


Johnson, 
et al. 
(1980) 

Arrow 

71.2 - 

43 

1.65 

0 T.E. Flap 

Manro 

(1983) 

Cropped 

Arrow 

71.2 - 

43 

1.65 

0.1 Twist 

Twist/ Camber 


8.3 


“ 

Thin 

- 

Full 

11.9 

0 

0.4 




Thin 


Full 

16 

0 

0.4 


Harrison Cranked 70 58 48 2.24 0 L. E. Droop 0-15 58 - Thin - Part 20-25 0 0.65 

(1982b) Arrow + 


+ AFT I /F- 111 Planform and Camber 
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TABLE 1.- CONTINUED 


(g) all planforms, conical camber 


Planform 





Conical Camber 



Thickness 

Vortex Flow 



Reference Class 

a 2 

He 

AR 

X 

Class 

A hl 


Wing Fusl 

Span a 

p 

M 


Johnson, 
et al. 
(1980) 

Delta 

74 


0 

1.15 

0 

Wentz 




Thin 


Full 

10-31 

0 

0 

Kuhlman 

(1978) 

Delta 

74 

- 

0 

1.15 

0 

Wentz 



- 

Thin 

- 

Full 

10-31 

0 

0 

Lamar 

Delta 

74 

- 

0 

1.15 

0 

Wentz 

- 

- 

- 

Thin 

- 

Full 

15-25 

0 

0 

(1978) 


80 

" 

" 

0.71 

" 

Barsby, p=0.2 

" 

” 

" 

" 

" 

" 

10-25 

" 

*' 








Bars by, p=0.4 







15-25 



Luckring, 

Delta 

76 

- 

0 

1.00 

0 

Barsby, p=0.4 

- 

_ 

_ 

Thin 

_ 

Full 

14 

0 

0 

et al. 
(1985b) 







Barsby, p=0.6 







20 



Reddy 

Delta 

76 

- 

0 

1.00 

0 

Nangia Wing E 

- 

- 

- 

Thin 

- 

Full 

10-30 

0 

0 

(1979) 



" 

" 

" 


Squire Wing 1 

*' 

" 



" 

" 

10-40 

" 







•• 


Squire Wing 2 
Squire Wing 3 







15-40 







•• 


•• 

Squire Uing 4 
Squire Wing 5 







30-40 



Reddy 

Delta 

76 

- 

0 

1.00 

0 

Nangia Wing B 

- 

- 

- 

Thin 

- 

Full 

24,30 

0 

0 

(1981) 


" 

*’ 

” 

" 

" 

Squire Wing 4 


” 


*’ 


" 

15-40 

” 

" 




“ 

" 

" 

" 

Squire Wing 5 



•* 

" 



20,30 

" 

*' 








Squire Wing 6 







20,30 



Tinoco, 

Delta 

74 

- 

0 

1.15 

0 

Linear Twist 

10 

- 

_ 

Thin 

- 

Full 

16-30 

0 

0 

et al. 




*• 

" 

“ 

■* 

20 






17-30 


" 

(1979) 







Wentz/Tabbed 







15-30 
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TABLE 1.- CONCLUDED 
(h) all pi an forms, vortex- flap camber 


Planform 





Vortex Flap Camber 


Thickness 

Vortex Flow 


Reference Class 

a 2 

A te 

AR 

\ 

Class 

A M 6 2 

Wing Fusl 

Span a p 

M 


Erickson 

Delta 

65 

- 

0 

1.87 

0 

Conical 

-15-40 

74 


- 

Thin 

- 

Full 

15 

0 

0.6 

(1983) 


" 

" 

*' 

" 

" 

i. . 

15 

" 



" 

" 

" 

11 

" 

" 



" 

" 


" 

" 

" 

20 

" 


" 


" 

" 

15-20 


" 



" 



" 


“ 

30 





" 

" 

15 


0 



" 

'* 


" 


Inverse Taper 

15 

” 




" 

Part 

" 


0.6 


Cranked 

72 

62 

" 

1.68 

" 


15 

77/70 

•* 

" 


Full 

10,15 

" 

0.6 


Delta 

















Erickson 

Delta 

55 

_ 

0 

2.80 

0 

Tapered Chord 

0-30 

53 


0-10 

Thin 

- 

Full 

10-30 

0-20 

0.4 

(1985) 


















Frink 

Delta 

74 

_ 

0 

1.15 

0 

Constant Chord 

10 

74 


- 

Thin 

- 

Full 

14 

0 

0 

(1982) 



" 

" 

" 

" 


” 

" 



" 

" 

Part 


" 

" 



" 

" 

" 

" 

“ 

" 

0-20 

" 


" 


" 

" 



0.3 




" 

" 

" 

'* 

Conical 

20 

79 




" 

Full 

11-20 



Frink 

Gothic 

60 

74 

0 

1.36 

0 

Gothic 

30 

74 


_ 

Thin 

Thin 

Part 

11-19 

0 

0 

(1985) 




" 

" 

" 

” 

40 




" 

" 


15-21 

" 

" 




■ 

" 

" 

" 


30 





Thick 

Full 

14,15 



Harrison 

Cranked 

70 

58 

48 

2.24 

0 

Tapered Flap A 

15-45 

58 


- 

Thin 

- 

Part 

15-25 

0 

0.65 

(1982) 

Arrow + 






Tapered Flap B 

25-35 

” 



" 








** 

" 

" 

” 

" 

Tapered Flap C 

30 






*■ 







" 


" 


Tapered Flap D 

30 

" 



" 









" 


" 


Tapered Flap E 

30 

" 














" 


Twisted Flap B 

35/25 

" 






20,25 





" 

" 

" 

•• 

" 

" 

25/35 

" 


" 


" 


" 











40/30 





" 







" 





" 

45/30 







" 

" 







** 



35/45 

" 













" 



" 

45/35 





" 


20 

" 





" 




Apex 2 , Flap A 

30 







20-25 








" 


Apex 3, Flap A 

'* 







. 






65 

55 

1.64 


Tapered Flap A 

** 





" 


20,25 





" 

" 




Tapered Flap B 














58 

48 

2.24 



0-30 




" 



25 



Harrison 

Cranked 

70 

58 

48 

2.24 

0 

Flap D, 

35 

58 


0-23 

Thin 

- 

Part 

20-25 

0 

0.65 

(1982b) 

Arrow + 





" 

L. E. Droop 











Hoffler 

Delta 

65 

_ 

0 

1.87 

0 

Conical/Tabbed 

0-30 

74 


0-30 

Thin 

- 

Full 

16 

0 

0.1 

(1985) 


















Luckring, 

Delta 

72 

- 

0 

1.30 

0 

Conical 

13 

78 


- 

Thin 

- 

Full 

14 

0 

0.2 

et al. 


74 

" 

" 

1.15 

" 

" 

0-25 

79 


" 

*■ 


" 



0 

(1982) 




” 

" 

" 


20 

76 - 

-80 


*' 

" 

" 

" 


" 



*' 


" 

" 

" 

" 

" 

79 



" 


" 

11-20 










Constant Chord 


74 



" 


" 

14 

" 

0.3 

Reddy 

Delta 

74 

- 

0 

1.15 

0 

Constant Chord 

-30 

74 


- 

Thin 

- 

Full 

15-20 

0 

0 

(1981c) 


" 

" 

" 

" 

" 

Tapered Chord 

24 

71 

4 

” 

" 

" 

" 

15-35 

" 

" 









-60-35 







20.6 
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A = 55° 85 c 


A = 0.25 -0.4 


\ 

a= 2° -50° 

\ 


, P = 0° -15° 

\ - 

A 

1 M = 0-0.99 

0 


M : 


Arrow, Diamond 
A te = 55 —51 



a =7° -40° 
P= 0° 

M = 0 — 0.76 


Cranked wings 



A = 0, 0.2 


(a) flat wings 


Conventional 



Twist 

Twist /camber 
Conical camber 
T.E. Flap 
Smooth onflow 


Vortex Flaps 



Conical 
Inverse taper 
Constant chord 
Swept, tapered 
Tabbed 


JCjO 35© 

0 

0 - 0.6 


= 5° -30° 
= 0 

= 0 , 0.6 


= 10-40 
= 0 

= 0 - 0.6 


(b) cambered wings 

Figure 1.- Converged FVS parameter space. Published results 



Rectangular 
AR = 0.2 - 12 



Thick arrow 
t/c — 0.03 


|3=0 
M= 0.4 


Raked side edge 
5 + = -3-13 



^ V 


a= 23° 

P - 0 
M = 0 


Thick wing and body 
t/c -0.03 


a= 16° 

p=o 

M = 0.4 


(c) other wings 
Figure 1.- Concluded. 



Vortex modeling - Real flow effects 

• 3D pressure field 
• Accurate forces and moments 
• High -a vortex lift loss 

• Vortex curvature 

• Vortex crowding 

• Viscous flow extensions 

• Vortex breakdown 

• Secondary separation 


Figure 2.- General solution features. 
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Wing 

(b) vortex flow 
Figure 4.- Concluded. 



Near wake 


Suction Analogy 















Available Solution Sheet Transfer Sheet Transfer 

(Partial restart) and Rotation 

Figure 10.- Sheet rotation concept with partial restart. 


Starting Solution 



Converged 

Solution 



It 


A6f = 5° (small) 



(a) rotation; a = 65°, A hl = 74% « = 16° s M = 0 


Figure 11.- Effect of vortex sheet manipulation. 
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Vortex sheet 



(b) scaling; a = 74°, a = 21°, M = 0 
Figure 11. Concluded. 


Desired FVS solution ; 
74° delta w/fuseiage 




Figure 12.- Application of expanded restart to complex geometries. 



Starting jV^ 'P'/ V °° | Starting 



Figure 13.- Multiple vortex sheet partial restart. 


• Attached flow analysis 

• a so = a for smooth onflow 

• Seek solution for small Aa 

• Aa = a-a so 

• Short near wake 

• Vortex starting solution 

• size by conical flow at Aa 

• rotate by leading edge deflection 



Barsby (1974) 


Figure 14.- A direct approach for large camber. 
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0 0.2 0.4 0.6 0.8 1.0 1.2 

y/s 

(a) p = 0.4, a ~ 14°, a so « 9° 



0 0.2 0.4 0.6 0.8 1.0 1.2 

y/s 

(b) p = 0.6, a ~ 20°, « so ~ 14° 


Figure 15.- Converged results at trailing edge. AR = 1. 
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